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The discovery of cosmic acceleration has presented a unique challenge for cosmologists. As ob- 
servational cosmology forges ahead, theorists have struggled to make sense of a standard model 
that requires extreme fine tuning. This challenge is known as the cosmological constant problem. 
The theory of gravitational aether is an alternative to general relativity that does not suffer from 
this fine-tuning problem, as it decouples the quantum field theory vacuum from geometry, while 
remaining consistent with other tests of gravity. In this paper, we study static black hole solutions 
in this theory and show that it manifests a UV-IR coupling: Aether couples the spacetime metric 
close to the black hole horizon, to metric at infinity. We then show that using the Trans- Planckian 
ansatz (as a quantum gravity effect) close to the black hole horizon, leads to an accelerating cosmo- 
logical solution, far from the horizon. Interestingly, this acceleration matches current observations 
for stellar mass black holes. Based on our current understanding of the black hole accretion history 
in the Universe, we then make a prediction for how the effective dark energy density should evolve 
with redshift, which can be tested with future dark energy probes. 



I. INTRODUCTION 

The discovery of recent acceleration of cosmic expan- 
sion was one of the most surprising findings in modern 
cosmology 0, . The standard cosmological model (also 
known as the concordance model) drives this expansion 
with a cosmological constant (CC). While the CC is con- 
sistent with (nearly) all current cosmological observations 
psf , it requires an extreme fine-tuning of more than 60 
orders of magnitude, known as the cosmological constant 
problem 

In the context of the concordance cosmological model, 
there are (at least) three different aspects of the CC prob- 
lem. For decades, physicists worried about why the value 
of the cosmological constant/vacuum energy seemed to 
be nearly zero by particle physics standards (known as 
the old CC problem)^, and the conventional wisdom was 
that it should vanish exactly, as a result of a yet unknown 
symmetry of nature. The accelerated cosmic expansion 
has thus challenged us to address this question on two 
new fronts. First is the new CC problem: why is the vac- 
uum energy density so close to zero, but non- vanishing? 
Second is the coincidence problem: Why did the dark en- 
ergy dominance and structure formation happen at ap- 
proximately coincident times? 

The race is on to simultaneously address these three 
questions. A popular alternative approach to the cos- 
mological constant is a model that modifies Einstein's 
theory of gravity. Traditionally, this involves adding 
higher order curvature terms to the geometric side of Ein- 
stein's equation. However, in ||^, one of us proposed a 
novel approach to modified gravity. This model intro- 
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duces gravitational aether, as a necessary ingredient to 
decouple the quantum field theory vacuum from gravity 
while simultaneously satisfying other tests of gravity. Un- 
like many models of modified gravity, the gravitational 
aether model modifies the energy-momentum content of 
the spacetime, instead of adding higher order curvature 
terms. 

In this model, the right hand side of the Einstein field 
equation is modified as: 

{SnG'r'G^, = T^, - -T^g^, + p\u'^K + 9^^), (1) 

where G" is 4/3 times the Newton's constant, and p' and 
u'^ are aether pressure and four-velocity that are fixed by 
requiring the conservation of the energy-momentum ten- 
sor, T^i,, and the Bianchi identity. As argued in |^, while 
consistent with precision tests of gravity, this theory is 
preferred to the standard model by the combination of 
cosmological observations (with the notable exception of 
■*He primordial abundance). 

In this paper, we pursue a detailed understanding of 
static spherical black hole solutions in the gravitational 
aether theory. The solution we find is, at first glance, a 
perturbed Schwarzschild metric. However, upon closer 
inspection we find that this perturbation is divergent 
both near to and far away from the horizon (where we 
refer to an infinite redshift surface as a horizon). Thus 
the static solution in the presence of gravitational aether 
is fundamentally different from Schwarzschild, which can 
be characterized as a UV-IR connection: the metric near 
and far from the horizon is set by the same integration 
constant. Here, we will explore possible meanings of this 
property, and whether the cosmological behavior is set by 
a Trans- Planckian ansatz close to the black hole horizon. 

We note that the static black hole solution found here 
also applies to the cuscuton models which have 
the same energy-momentum tensor as the gravitational 
aether in the limit of vanishing cuscuton potential. 
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In Sec. H, we introduce our gravitational aether black 
hole solution. We describe the properties of the solution, 
including a preferred coordinate system and the location 
of the event horizon. We also establish asymptotic prop- 
erties of the black hole, which are characterized by the 
same integration constant both close in and far away from 
the horizon of the black hole. 

Sec. Ill explores the Trans-Planckian ansatz, as a way 



to fix the aforementioned integration constant, through 
quantum gravity effects close to the horizon. We sug- 
gest a way to connect the presence of black holes to the 
existence of a pervasive pressure that behaves like dark 
energy on cosmological scales. 
Sec. IV, we present 
many such black 
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tion that 



a study of the contribu- 
holes would make to the 



global/cosmological structure of space-time, while Sec. 
|y| provides a census of average black hole mass through 
cosmic history, which translates into a prediction for the 
history of cosmic acceleration. 

Finally, we will discuss open questions and future 
prospects in Sec. VI. 

Throughout the paper, we use the natural Planck 
units: h ~ c = Gn ~ ks — 1- Moreover, we will re- 
place pressure p' by 3p/4 in Eq. (|^), so that the vacuum 
field equations for the aether theory resembles general 
relativity sourced by a perfect fluid with pressure p and 
zero density. 



II. BLACK HOLE IN GRAVITATIONAL 
AETHER 

We find a solution for the static black hole in the Gravi- 
tational Aether model using assumptions similar to those 
that lead to the Schwarzschild solution. Namely, we as- 
sume a spacetime with no matter content, and we assume 
spherical symmetry. Given that the aether takes fluid 
form, the metric in this model is the same as the general 
static, spherically symmetric metric that describes the 
interior of a star, as modeled by a perfect fluid. The only 
notable divergence from the star model is the absence of 
a matter density, leaving an energy-momentum tensor of 
the following form: 

Tf^iy = p{u^,u,y + gi_,^). (2) 

We find the following metric 

ds^ = -e'Ue+l^^'^^ \r^+r'dn^. (3) 

With components obeying the following differential equa- 
tion, known as the Tolman-Oppenhcimer-Volkoff equa- 
tions l|§l: 
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FIG. 1: Function f{r) [Eq. @] as a function of the distance 
from the Schwarzschild radius {— 2m). The deviation from 
the Schwarzschild metric is proportional to pof(r), where po is 
the integration constant. If po is small, as we argue in Section 
the corrections only become important at the horizon, and 
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on cosmological scales. 



We see immediately that exp(0) and p are inversely re- 
lated: 



P = Poe ^, 



(6) 



where po is an integration constant. Notice that Eq. (|^) 
is equivalent to the condition of hydrostatic equilibrium 
for aether, and is valid independent of the assumption of 
spherical symmetry, for any static spacetime ||29|| . Now, 
we may rewrite Eqn. 0: 
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(7) 



We can solve this equation by noting that it is a first- 
order inhomogeneous linear differential equation in e*^, 
with the standard solution: 
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To put this into more familiar terms, we can set the 
constant, so that we recover the Schwarzschild solution 
as Po — > 0: 



where /(r) is given by: 
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and is shown in Fig. (|^). In the hmit where r is large 
{r :S> m): 



fir) = Y + ^rnr + 0[m'^]. 
while close to the "Schwarzschild horizon" we find: 



(11) 



fir) = -i 



\/—2m + r 



+ 0[m^/^ir^2m)^/^]. (12) 



Thus the correction to the Schwarzschild metric domi- 
nates in both UV and IR regimes (corresponding to close 
to and far from the BH horizon). This a nice tie, even 
for arbitrarily small values of the integration constant po. 
Therefore, a very suggestive conclusion is that, unlike in 
general relativity, the gravitational aether ties the forma- 
tion of black hole horizons to cosmological dynamics. 

But then, is there really an event horizon for this space- 
time? Looking at the trace of the Einstein's equation, we 
find that the Ricci scalar is proportional to the pressure 
of aether, p, which is in turn inversely proportional to 
the 00 component of the metric, e"^. We define the sur- 
face where e*^ — > as the black hole horizon. Therefore 
the pressure at the horizon, and thus the Ricci scalar, 
goes to infinity (p oc TZ oo) implying that this surface 
coincides with a real metric singularity (as opposed to a 
coordinate singularity). 

It appears that we have established that any static 
event horizon in a theory of gravitational aether (like the 
one we have modeled) coincides with a real metric singu- 
larity. In a traditional formulation of general relativity, 
such a scenario may be given to ambiguous physical inter- 
pretation. Cognizant of the fact that a modified gravity 
will display properties divergent from traditional relativ- 
ity, we expect that such a picture is best contextualized 
by a more comprehensive theory of quantum gravity. 

Indeed, any process (for example, quantum gravity) 
that alleviates/regulates metric singularities, will in- 
evitably remove event horizons from the theory of grav- 
itational aether. In other words, static event horizons 
cannot exist in a UV completion of gravitational aether. 
This is independent of the assumption of spherical sym- 
metry, and only relies on the aether hydrostatic equilib- 
rium condition (^. However, we note that, as the sin- 
gularity is a null surface, the spacetime does not violate 
the weak cosmic censorship principle. 

Back to the spherical aether black hole spacetime (||), 
we now notice that the static metric solution is only well- 
defined for r > 2m, as the solution becomes complex 
inside the Schwarzschild radius, r < 2m. More surpris- 
ingly, for negative values of po, unlike a Schwarzschild 
black hole, a free-falling observer can reach this bound- 
ary within a finite coordinate time. The reason is that 
the redshift of a static source at the Schwarzschild radius 
is now finite as seen by distant observers [pOl : 
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As to what happens inside r = 2m, it is clear that our 
current choice of coordinates do not give us a physical 
metric for r < 2m. However, is it possible that with 
an appropriate choice of coordinate, we can analytically 
continue the static solution beyond the Schwarzschild ra- 
dius? Indeed, we can define a new radial coordinate: 



A = / dr'^/g^ = / i 

J 2m J 2m 



-1/2 



= 2 [2m,(r - 2m,)]^/^ + O [(r - 2m)3/2„j-i/2l ^^-^^^ 

which is equivalent to the constant-time proper radial 
distance. In terms of A, the metric takes the form: 



-e^*dt^ + d\^ +r{\)^dn^ 



(15) 



where 



e't' = -327rpom2 -f — + 0\pq\^ , A^m-^/^l, (16) 
Am 

\2 

r(A) = 2m, + + 0\\^/m\ (17) 

8 m 

In other words, the metric is analytic and real in terms 
of the new radial coordinate. A, at and beyond the 
Schwarzschild radius, which now corresponds to A = 0. 
Moreover, a static event horizon, which as we argued cor- 
responds to a real curvature singularity, now exists for all 
(small) values of po, as e"^ = at: 



(18) 



In the next section, we study the implications of this 
solution for cosmology. However, we shall postpone a full 
investigation of the causal structure of this spacetime, 
as well as its possible analytic continuations, to future 
studies. 



III. TRANS-PLANCKIAN ANSATZ AND 
COSMIC ACCELERATION 

In the last section, we saw that within spherical space- 
times in the gravitational aether theory, the integration 
constant po ties the geometry close to the horizon to the 
geometry at infinity. While, in the classical theory, p^ is 
an arbitrary integration constant, here we speculate that 
its value is fixed by quantum gravity effects, especially 
since the horizon is now a curvature singularity, where 
quantum gravity effects should become important. 

We first note that the temperature of sources that fall 
through the Schwarzschild horizon, as seen by distant 
observers, approaches the Hawking temperature 0: 
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FIG. 2: Predicted large distance deviation from the vacuum 
Schwarzschild solution for 1, 10, and 100 Mq black holes, 
based on the Trans-Planckian ansatz. Here, we assumed 6p = 
100 in Eq. (^2|) for non-rotating black holes to find po, which 
is then plugged into Eq. to find the metric. As pointed 
out in the text, the corrections become important on today's 
cosmological horizon scale for solar/stellar-mass black holes. 



Furthermore, we assume that the maximum rest-frame 
temperature of sources is comparable to the Planck tem- 
perature (or one in Planck units): 



0[l]. 



(20) 



Here, 9p is a dimcnsionless constant that measures Tmax 
in units of Planck temperature, which we shall call the 
Trans-Planckian parameter. We then adopt the Trans- 
Planckian ansatz, which is the idea that the maximum 
redshift at Schwarzschild radius (Eq. |l^) is roughly set 
by the ratio of the Planck to Hawking temperatures: 



1 + 
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With this ansatz, we further sec that: 
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= o[i], 



(21) 



(22) 



(23) 



i.e. the event horizon is roughly a Planck length away 
from the Schwarzschild radius. Equivalently, the short- 
distance aether corrections to the Schwarzschild metric 
only become important at about a planck distance from 
the horizon/singularity, which is a reasonable expecta- 
tion from a possible quantum gravitational mechanism. 

While this may imply that tests of strong gravity close 
to the horizon of a black hole may have a hard time 
testing the influence of aether on the spacetime metric, 
the Trans-Planckian ansatz has a curious prediction for 
the numerical value oipo, i.e. aether pressure far from the 



black hole. Comparing the scale of po with the density 
(~ — pressure) of the cosmological dark energy, p\: 
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(24) 



where we assumed Q\ = 0.7 and Hq = 70 km/s/Mpc. 
The resulting deviation from the Schwarzschild metric is 
shown in Fig. (|^) for stellar mass black holes. 

This leads us to a very interesting possibility, which 
was first conjectured in H: that the formation of stellar- 
mass black holes could trigger the onset of cosmic ac- 
celeration, especially since aether and dark energy have 
similar pressures, assuming that the aether pressure is 
set by the Trans-Planckian ansatz for stellar mass black 
holes. To see this, we can explicitly compare the black 
hole spacetime (Eqs. ^ and ^T|) far from the black hole 
(r > m): 



ds^ = -(1 + 2TTpar'^fdt^ + dr^ 
with the de-Sitter spacetime: 



(25) 



ds^ = -{l-STrpfj-^ /i)dt^ + {lS-KpKr^- /iy^dr'^Wd^f. 

(26) 

We thus notice that non-rclativistic particles close to 
the origin, but far from the black hole horizon (2m <i; 
f ^ bo|~^^^) see the same Newtonian potential (or gu) 
in both spacetimes, if po = — 2pa/3. In other words, 
close-by non-relativistic test particles (such as galaxies, 
stars, or other black holes) accelerate away from the ori- 
gin/black hole, similar to a de-Sitter space. Moreover, 
this acceleration will correspond to the current cosmolog- 
ical observations, if the mass of the black holes is roughly: 
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-1/3 



(27) 



So far, our solution has neglected the effects of black 
hole spin. Indeed, spin is expected in realistic black holes, 
which are fed by astrophysical accretion disks. For ex- 
ample, the dimensionless spin parameter, a* = a/m was 
recently measured for two stcllar-mass black holes, to be 
within the range 0.65 — 0.85 |0. In order to include 
this effect, we conjecture that po scales as Tfj (as sug- 
gested in 1^), for general black hole spins. This is justi- 
fied, as the Trans-Planckian ansatz is controlled by the 
Hawking temperature, Ty, while /(r) also depends on 
the surface gravity close to the black hole horizon, which 
is also proportional to Th- With this assumption, the 
scale-dependence should go as: 



Po{m,a^) 
Po(to,0) 



8 



-1/2 



(28) 



which is in the range 0.2 — 0.6 for a* = 0.65 — 0.85. 

While this paper only deals with static vacuum solu- 
tions, it was shown in jsj that for non-relativistic fluids 
(e.g. stars, planets): p' « — T^/4-|-const., i.e. the local 
matter density sets the aether pressure up to a constant. 
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One expects that the constant term would be set by the 
boundary conditions at infinity, or by cosmology. Al- 
ternatively, what we suggest in this section is that the 
boundary condition can be set at the horizons of the 
black holes. The fact that this can naturally explain the 
onset of cosmic acceleration is certainly very suggestive, 
but the best way to test this hypothesis is to see how/if 
this boundary condition can emerge from the process of 
(classical or quantum) gravitational collapse into a black 
hole. We leave this question to future studies. 

A further implication of this hypothesis is that so- 
lar/stellar mass is the minimum mass of black holes al- 
lowed in the model. A discovery of significantly sub- 
solar mass black holes (e.g. primordial black holes 
with M ^ Mq) could potentially rule out the Trans- 
Planckian ansatz, as it would imply much larger than 
observed cosmic acceleration rates for 6p ^ 1. 

Of course, we also need to patch together and coarse- 
grain individual black hole spacetimes into a de-Sitter 
space, in order to rigorously prove this correspondence. 
However, the above argument is already very sugges- 
tive, as long as there are many black holes within the 
cosmological/de- Sitter horizon, so that one can trust the 
above Newtonian argument. In the next section, we pro- 
vide an approximation to the cosmological spacetimc of 
multiple black holes. 



IV. GLOBAL CONTRIBUTION OF MULTIPLE 
BLACK HOLES 

In this section, we will seek an approach to approxi- 
mately find the spacetime of multiple black holes with 
gravitational aether, which can be used to describe an 
approximate FRW cosmology. Here, for simplicity, we 
focus on the quasi-static Newtonian regime, where we 
could assume hydrostatic equilibrium for aether in the 
vacuum (|^). For simplicity, we ignore the matter in- 
between black holes ||3l|, and assume that black holes 
are much farther apart than their horizon sizes, but are 
much closer than the cosmological horizon. In this limit, 
using Eq. ^ we have: 



Green's function, which can be found using the method 
of images (e.g. 



Gi5(x,x') 



1 



|x-x'| a;'|x-^x' 



(30) 



For n spheres (black holes) at positions x^ and with radii 
ai (= 2mi), we may expand this Green's function, up to 
first image, as 



G_d(x,x') = 



1 



= 1 |X' - Xj||X - Xj 



|x'— Xi I' 



■(x'-x,) 



O 



|Ax|3_ 
(31) 



which is a good approximation, as long as the distance 
between the spheres/black holes is much larger than their 
sizes. Now, using Green's theorem, we can find aether 
pressure in-between the black holes, in terms of the pres- 
sure on the surfaces of the spheres, pi's: 

Inp(x) - Inp = -— ^ * ds' • -^[lnp,(x') - Inp], 

(32) 

where Inp is the log of pressure at infinity, and §g ds' 
are surface integrals over the horizons of the black holes 
(assuming a flat geometry), while pi cx are fixed by 
the masses of the blacks holes, using the Trans-Planckian 
ansatz (p^). Since the Green's function (^Tj) is analo- 
gous to superposition of electrostatic potentials of point 
charges, we can use Gauss's theorem to evaluate the sur- 
face integrals: 



Inp(x) — Inp = 



ai[lnpi(x') - lap] 



(33) 



Now, using the assumption of statistical homogeneity, 
we expect the spatial/ensemble average of \i\p to be the 
same as In p. If we take ensemble averages of both sides 
of Eq. (33), this yields: 



VMnp = -V^t/) = 0, 



(29) 



\np 



{at In Pi 



(34) 



where the assumption of = is the equivalent of the 
Poisson equation in Newtonian gravity, for zero matter 
density (which also applies to aether). We thus see that 
fixing the aether pressure in the vicinity of black holes, 
through the Trans-Planckian ansatz (p^, is equivalent to 
solving the Laplace equation (^9|) with Dirichlet bound- 
ary conditions at (or close) to the horizon of the black 
holes [|2|. 

This problem is analogous to finding the electrostatic 
potential between multiple conducting spheres, which 
can be solved using the Green's function for the appropri- 
ate geometry. For a single sphere of radius at the origin 
(and in a flat space), there is an exact expression for the 



or alternatively: 



256Tr^0pm3 



In TO* 



{mi In TOi) 



(35) 

where we used pi cc m^ '^ and Oi = 2mi, as well as Eq. 
(^). In other words, in the presence of multiple black 
holes, the mean aether pressure, and thus FRW cosmol- 
ogy, is set by m*, which is the mass- weighted geomet- 
ric mean of black hole masses. Subsequently, the cor- 
respondence of this mean aether pressure with an effec- 
tive Dark Energy or cosmological constant density was 
demonstrated in the last section. 
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Furthermore, taking the Laplacian of Eq.(^), we can 
find an equation for the perturbation of effective Dark 
Energy, for sub-horizon perturbations (but on scales 
larger than the size of the blacks holes): 



(36) 



where feE.off is the overdensity of the effective dark en- 
ergy, while pbh is the black hole density. Eq. (|3^) can be, 
in principle, used to track cosmological structure forma- 
tion and the impact on CMB anisotropies (through the 
Integrated Sachs- Wolfe effect), but we postpone a study 
of these effects to future work. 

In the next section, we will provide a quantitative pic- 
ture of how the cosmic history of accretion into stellar 
and super-massive black holes (or active galactic nuclei) 
leads to an estimate of function of redshift. and 

its implications for the effective dark energy scenarios. 



V. COSMIC HISTORY OF BLACK HOLES AND 
COSMIC ACCELERATION 

An up-to-date inventory of cosmic energy at the 
present day, including the contribution from stellar-mass 
and supermassive black holes, is provided by [p^ . In or- 
der to measure m* (z) we need to take this a step farther 
and understand the mass distribution of such black holes, 
and their redshift evolution. 

The mass distribution of stellar-mass black holes is 
not well-determined observationally, but estimates are 
that it is fairly broad, with a mean of around ~ 7Mq 
p3| , . We base our calculations on the theoretical pre- 
dictions of , which show that the distribution can be 
approximately represented by a power-law such that the 
number density of black holes decreases by a factor 5 
between M = 3Mq and M = 15Mq. Assuming this 
distribution, the average black hole mass is 8.2Mq. We 
will treat the uncertainty in this distribution by varying 
the slope sufficiently to alter this mean mass by zLIMq. 
To determine the redshift evolution in black hole abun- 
dance, we use observations of the cosmic star formation 
history, from There is significant uncertainty in the 
shape of this history; however it must obey the integral 
constraint that the total stellar mass density today be 

= 0.0027pcrit = 3.67 X IO^A/qMpc^^^ which is known 
to a precision of ~ 30 per cent We will therefore 

normalize the black hole number density at z = to 



1.46 X 10*^ Mpc"^ at z = [0. We assume that changes 
to the initial mass function do not significantly alter the 
shape of the star formation rate density evolution, but 
primarily affect the number of black holes formed. By 
default we assume a Kroupa IMF 1 1 7 1 , which is the "sec- 
ond model" considered in For this choice, 0.19 per 
cent of stars formed end up as black holes; a more use- 
ful number is that for every solar mass of stars formed 
0.0025 black holes are created. These numbers change by 
less than 5 per cent if we assume a Chabrier IMF p8| ; 
we expect therefore the uncertainty on the normalization 



of the black hole mass function to be dominated by the 
30 per cent uncertainty in the present day stellar mass 
function. Note, however, that a pure Salpeter IMF |l^ 
would produce significantly fewer black holes, only 0.0013 
for every solar mass formed. 

We base our estimate of the supermassive black hole 
mass distribution on observations of the quasar lumi- 
nosity function. This requires assumptions about the 
lifetime and obscuring column density of quasars; for 
this we adopt the model of |Q who describe a merger- 
driven scenario of black hole growth. Using this model, 
the z = mass density of supermassive black holes is 
2.9tl'^ X lO^MgMpc"^. This is somewhat smaller than 
the value of 5.4 x lO^MoMpc"^ determined from the cor- 
relation between black hole mass and bulge luminosity 
pl| , |2^ , as computed by ||l^. However, the uncertainty 
on the latter is a factor of two, and a lower value of 
3.4 X IO^MqMpc"^ is obtained if one uses the cor- 
relation with velocity dispersion for early type galaxies 
p3| , rather than luminosity. 

With this in hand we are able to compute the expected 
m*(z), and this is shown in the bottom panel of Figure]^. 
Our best estimate of the local, mass-weighted geometric 
mean of black hole masses is m*(0) = 12.7M0. The 
dashed lines represent the range of uncertainty on this 
z = normalization. A larger value of m» is obtained 
by reducing the contribution of stellar-mass black holes 
(assuming the local density is 30% lower than our fidu- 
cial model, and assuming the mass distribution is more 
steeply weighted to lower masses, so the average mass is 
7.2M0), and increasing the contribution of supermassive 
black holes (by increasing the z = space density within 
the Icr uncertainty, to 5.2 x lO^MgMpc-^). This yields 
TO*(0) = 24.7M0. Pushing the numbers in the opposite 
direction, we obtain m*(0) = 1O.5M0. Using Eq. ( [27| ) 
for the current effective density of dark energy, and ignor- 
ing the spin of black holes, this range in to*(0) translates 
to a range for the Trans-Planckian parameter 9p: 



= (0.4-5) x 10^ 



(37) 



We can consider spinning black holes, using our scaling 
argument from Section 111 and taking a nominal value 
of a* = 0.75. This implies a lower range for the Trans- 
Planckian parameter, 9p = 20 — 300, in order to match 
the current rate of cosmic acceleration. The fact that 
6p ^ 1, further justifies a Trans-Planckian, or quantum 
gravitational origin for the observed "dark energy phe- 
nomenon" . 

The evolution of the stellar-mass black hole mass den- 
sity is dependent upon the shape of the star-formation- 
rate density plot from ||l^. To consider the effect of this, 
we construct two star formation histories that are consis- 
tent with those data within the Icr error bars, but which 
produce as many stars as possible at either high redshift 
(z > 1) or at low redshift (z < 0). We still renormal- 
ize this to match the local stellar mass density. These 
extremes are shown in Figure ^ as dashed lines. The 
evolution of the supermassive black hole distribution is 
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We can define a mean equation of state as: 




Redshift 

FIG. 3: Bottom panel: The mass-weighted geometric mean 
of black hole masses, m*, in units of AIq as a function of red- 
shift. Our fiducial model (solid, black line) assumes our best 
estimates of the mass distribution evolution of the black hole 
mass distribution. Dashed lines indicate the range of uncer- 
tainty expected due to the unknown relative contribution of 
supermassive and stellar-mass black holes at z = 0, while 
the dotted lines represent the uncertainty in the shape of the 
star formation density evolution from jl^]. Top panel: The 
prediction of the equation of state parameter ■w{< z) from 
Equation for the same models. The dashed area shows 
the region excluded at 68% confidence level for this parame- 
ter, as measured from independent observations |2a]. 



very model dependent, and not well constrained. We note 
that the two different predictions shown by ||2^, which 
make different assumptions about the quasar space den- 
sity evolution at z > 2, have a subdominant effect on the 
predictions shown here, relative to the other uncertain- 
ties considered. 

Within an effective dark energy description of FRW 
cosmology, a fixed dark energy equation of state, w, im- 
plies that dark energy density evolves as {1 + zY^^^^\ as 
a function of redshift, z. The effective equation of state 
(which is simply a way to parameterize cosmic expansion 
history) is observationally constrained to 



w{z) -1.06 ± 0.14 + (0.36 ± 0.62) 



1 



(38) 



at 68% confidence level, based on cosmic microwave back- 
ground, baryonic acoustic oscillations, and supernovae 
la observations, assuming a spatially flat cosmology [ p5[ , 
and a linear dependence of w{z) on the cosmological scale 
factor = (1 -I- z)~^. 



1 -I- w{< z) = 



_ 1 ln[pDE,eff(-z)/yODE,eff(0)] 



ln(l + z) 
ln[m,(z)/TO*(0)] 
ln(l + z) 



(39) 



since, /3DE.off(-z) esc m~'^(z), as we saw in the last sec- 
tion. We show this estimate of w for the models described 
above, in the top panel of Figure |^. Our fiducial model 
predicts a value of w that deviates from — 1 by less than 
5 per cent out to z '-^ 2, but predicts it should reach 
w = —0.8 by z = 3. There is considerable uncertainty 
on this, however, due both the unknown distribution of 
black hole masses at z = (dashed lines) and the un- 
known shape of the star formation rate density evolution 
(dotted lines). 

While most these models are consistent with the cur- 
rent bounds on the effective dark energy equation of state 
(using Eq. |3|): 

w{< z) = -1.06±0.14+(0.36±0.62) 1 - 

(1 + z) ln(l + z) 

(40) 

stage IV dark energy missions, as quantified by the dark 
energy task force report |^ , are expected to have percent 
level sensitivity to ?«(< 1 — 3), and thus should be able to 
distinguish the aether model with these m*(z) histories 
from a cosmological constant. 



VI. CONCLUSIONS AND FUTURE 
PROSPECTS 

We have shown that static black hole solutions exist in 
the gravitational aether model of Q . The model is an at- 
tractive alternative to the cosmological constant, which 
does not suffer from the tremendous fine-tuning problem 
of vacuum energy in standard model. We find that in 
the presence of a gravitational aether, the Schwarzschild 
black hole is sufficiently perturbed so as to result in 
a Trans-Planckian connection between physics near the 
black hole horizon and cosmology. This could be a phe- 
nomenological product of quantum gravity, and it nat- 
urally explains the present-day acceleration of cosmic 
expansion as a result of formation of stellar/solar-mass 
black holes. 

Indeed, the recent discovery of cosmic acceleration, or 
dark energy ^ might be the first concrete evidence for 
quantum gravity and/or Trans-Planckian physics. Fu- 
ture work may include an exploration of quantum prop- 
erties of this black hole solution. In particular, a natural 
next step would be to understand how quantum gravity 
can resolve the null singularity at the event horizon. 

As discussed in Sec. [II, another important question 



yet to be addressed is whether dynamical evolution could 
lead to the static solutions found in this work. While 
prior to formation of black holes, the integration constant 
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Po is set by large-scale conditions, as black hole horizons 
form, we speculate that the constant is instead set by 
conditions at the event horizon. In order to understand 
the causal transition between these two boundaries, and 
how fast the effect will propagate away from the black 
hole, a more complete dynamical picture is necessary. 

Furthermore, in the presence of multiple black holes 
with relative motion, the aether is expected to be locally 
dragged by different black hole horizons. However, for 
black holes at large separations compared to their hori- 
zon sizes and non-relativistic velocities (as expected in 
astrophysical situations), the perturbations to the static 
solution is expected to be small. 

To conclude, we would like to entertain the exciting 
possibility that the gravitational aether might provide 
a complete solution to the three aspects of the cosmolog- 
ical constant (CC) problem, as discussed in the Introduc- 
tion: 

1. Old CC problem: Gravitational aether theory de- 
couples quantum vacuum from geometry, which 
allows a nearly flat spacetime even in the pres- 
ence of large vacuum energy densities expected 
from the standard model of particle physics. The 
model makes specific predictions for physics at big 
bang nucleosynthesis and radiation-matter transi- 
tion era, which will be tested with precision cosmo- 



logical probes over the next decade ||5[] . 

2. New CC problem: Formation of black holes leads 
to a UV-IR coupling, which connects near-horizon 
Planck-scale physics to cosmology, and can natu- 
rally lead to cosmic acceleration, even in the ab- 
sence of a real dark energy component. 

3. Coincidence problem: As we showed in Sec. (^), 
the stellar mass black holes expected in standard 
star formation, can naturally lead to the observed 
present-day acceleration of the Universe. The com- 
petition between the contribution of stellar mass 
black holes, and super-massive black holes leads 
to an evolution of the effective dark energy den- 
sity, which can be tested with NASA's future Joint 
Dark Energy Mission ( JDEM) ||3^ or its European 
counterpart Euclid p% . 
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Here, distant observers are located at 



2m< r < (-po)"^''^- 

This is not a bad approximation since, as we argued in 
the last section, the effect of matter on the aether 
pressure is localized and does not extend into vacuum 
in the non-relativistic regime. 

Since pressure approaches po at several BH horizon 
radii for individual black holes, as long as the distance 
in-between black holes is much larger than their horizon 
radii, the exact radius at which the boundary condition 
is set is not important. 
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